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I. INTRODUCTION 

Many spatially extended systems that undergo a symmetry breaking instability close to a second-order critical point 
can be described by a real order parameter equation in the form of a Swift-Hohenberg type model, which has been 
derived in various field of nonlinear science such as hydrodynamics [1] , chemistry [2] , plant ecology [3] , and nonlinear 
optics [4]. 

A complex Swift-Hohenberg equation was deduced in the context of lasers [5-7] and optical parametric oscillators 
[ 8 ]. Moreover, to describe the nascent optical bistability with transversal effect in nonlinear optical cavities a real 
approximation has been deduced [9] from laser equations. This approximation allowed to predict stable localized 
structures and organized clusters of them [9] . A detailed derivation of this equation from first principles can be found 
in Ref. [5]. In this work, we show that this real modified Swift-Hohenberg equation (SHE) of the form 

dtu = r] + eu — u^ — ( 1 ) 

supports a curvature instability on localized structures that leads to an elliptical deformation that produces a rod-like 
structure. As the time evolution is further increases, the rod-like structure exhibits a transverse undulation and leads 
to the formation of invaginated structures. This structure is a labyrinthine pattern diefined by an interconnected 
structure where the field value is high. The outer part or complement to the invaginated structure corresponds to low 
held value. This behavior occurs far from any pattern forming instability and requires a bistable behavior between 
homogeneous steady states. In Eq. (1), u = u{x, y, t) is a real scalar held, x and y are spatial coordinates and t is time. 
The parameter y represents the external forcing held which brakes the rehection symmetry u —> —u. The bistability 
parameter is e. The coefficient ly in front of a diffusive term may change the sign and allows the pattern forming 
to take place [4, 10-13]. Depending on the context in which this equation is derived, the physical meaning of the 
held variable and parameters adopt different meanings, for instance, in cavity nonlinear optics u{x, y, t) corresponds 
to light held intensity, while parameters {y, e, v} are associated with the injection held, the deviations of the cavity 
held, and cooperativity, respectively ]9]. 

For certain range of parameter values, Eq. (1) exhibits stable circular localized structures, for y < 0 localized 
structures appear as isolated peaks of the held u{x,y,t), instead, for 77 > 0 localized structures are holes in the held. 
These localized structures have a hxed stable radius for each parameter value. Curvature instability of localized 
spot has been experimentally studied or theoretically predicted in magnetic huids [14], in liquid crystals [15, 16], in 
reaction-difusion systems [17-28], in plant ecology [29], in material science [30, 31], in granular huid systems and in 
frictional huids [32, 33], and nonlinear optics [34]. The hngering instability of planar fronts leading to the formation of 
labyrinth structures has been reported by Hagberg et al. [35]. In this manuscript we shall focus on circular localized 
states. 


II. STABILITY OF LOCALIZED SPOT 

Considering hxed parameter values, starting with an azimuthally symmetric localized structure. The structure is 
perturbed, this perturbation grows radially as shown in Fig. 2.1. The circular shape becomes unstable at some critical 
radius. The elliptical shape elongate into a rod like structures as shown in Fig. 1. This elongation proceeds until a 
critical size is reached beyond which a transversal instability onset the appearance of hngers near the mid section of 
the structure (see Fig. 2.3). The huger continues to elongate, and the amplitude of oscillation increases (Figs. 2.4 
and 2.5). The dynamic of the system does not saturate and for a long time evolution, the rod-like structure invades 
the whole space available in (x, 2 /)-plane as shown in Fig. 2.6. This invaginated structure is stationary solutions 
of the SHE. The dynamic described previously has been observed in cholesteric liquid crystals under the presence 
of an external electric held [15, 16], where an initially circular structure of cholesteric phase suffers from curvature 
instability, transversal oscillations and develops into an extended labyrinthine structure. The characterization of this 
dynamic is an open problem. 

For V = 2, the bifurcation diagram of the model Eq. (1) in the parameter space (e, y) is shown in Fig. 3. For 
e > 0 the system undergoes a bistable regime between homogeneous steady states. For e < 0, the system has only 
one homogeneous steady state. The curve Fi represent the coordinates of the limit points of the bistable curve 
{y± = ±2(e/3)^'^^). The threshold associated with a symmetry breaking instability is provided by the curve r 2 . 
The coordinates of the symmetry breaking instabilities thresholds are y± = ±\/(j^^ + 4e)/3(i/^ — 8e)/24. Curve Ta, 
built numerically separates the zone where localized structures are stable, //-zone, from the zone where they are 
unstable, /-zone. The transition from localized structures to labyrinthine pattern occurs when crossing from /-zone 
to //-zone through the Ta-curves indicated in Fig. 3. This transition occurs via fingering instability at the Ta-curves 
delimiting the parameter domain I and II. In the limit of the classical Swift-Hohenberg equation, 77 = 0, there is no 
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observation of fingering instability, instead the transition from II to I zones, localized structures only grow radially. 
As a result of boundary conditions destabilization of this structures into labyrinthine structures is observed, this is 
a size effect phenomenon. Contrary, for rj ^ 0 the transition from II to I zones of a localized spot is affected from 
curvature instability, giving rise to an unstable rod structure which exhibits transversal oscillations and develops into 
an extended labyrinthine structure. 

In what follows, we first study analytically the stability of a circular localized spot with respect to azimuthal 
perturbation. This mode analysis allows us to evaluate the threshold above which the transition from localized 
spot to a rod-like structure takes place. Then, we perform the linear stability analyzis of the rod-like structure and 
determine the conditions under which the transversal oscillations occur for the SH equation. 

Starting from a solution with rotational symmetry (i.e circular localized structure) u = Us{r — Ts) where r is 
the radial coordinate and characterizes the localized structure radius. Then we perturb the solution u{X, t) = 
Us{X) + W (r, TsjTq) with X is the relative position X = {r— (r^ +rQ{6, t))), where tq is the perturbed radius position, 
which accounts for the interface of the localized structure, 9 stands for the angular coordinate, and W 1 are 
corrections to the circular localized spot. Using polar representation of Eq. (1), considering the above perturbation 
and parameters in Eq. (1) at linear order in W one obtains 

CW = dtTodxUs + rj + eus - - lyV’^u^ - (2) 

where the lineal operator £ = —(e -|- 3Ug — — V"^). which is a self-adjoint. Assuming that the radius of the 

localized structure is sufficiently large, the operator £ which is explicitly dependent on the radial coordinate can be 
approximated by a homogeneous operator in this radial coordinate 


1 

• 

8A.U. 

2 

• 

3 

% 

4 

% 

5 

% 

6 

N 

7 

*-►x 

8 

\ 

9 

\ 


FIG. 1. Temporal evolution for; 1) t = 0; 2) t = 125; 3) t = 175; 4) t = 225; 5) t = 275; 6) t = 340; 7) 350; 8) 360; 9) t = 400, of 
a localized spot into an elliptical deformation and then to a rod-like structure for Eq. (1) with parameters: rj = —0.065; e = 2.45; 
1 / = 2.0. Minima are plain white. The image corresponds to a zoom of 16 x 16 points of a 512 x 512 point finite-difference 
simulation. 
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-e - iui + v[dl + —d; 



( 3 ) 


this operator possesses a neutral mode, i.e., zero eigenvalue with the eigenfuction dxUg. This approach allows us to 
perform analytical calculations, which are not accessible when the operator is inhomogeneous. Using this approach 
the solvability condition (see textbook [12] and references therein) yields 


6/3 


2 / 3 , 


dtro = -A^dgro + -^dgroidgroY - -^dgro + -^{dgro) 


( 4 ) 


where 


a {dxxUs\dxxUs) 

{dxUs\axUs) 

and 

A (i/ - 2/3). (6) 

Then the held ro(d,t) satishes a nonlinear diffusion equation. The hrst, the second, the third and the last term on 



FIG. 2. Transition from a single localized spot to invaginated pattern. Temporal evolution with Neumann boundary conditions 
and with the same parameters as in Fig. 1. 1) t = 0, Localized spot, 2) t — 600, rod-like structure, 3) t = 1900, transverse 
undulation of the rod-like structure 4) t = 2800, and 5) t = 3700, localized transient patterns, and 6) t > 15000, stationary 
invaginated labyrinth pattern. Minima are plain white and the mesh integration is 512 x 512. Simulation done with finite- 
difference method 




















5 


the right side of Eq. (4) account for the linear diffusion, the non-linear diffusion, the hyper-diffusion and nonlinear 
advection, respectively. When A < 0 the localized structure is stable, and for A > 0, the localized state is unstable 
as result of the curvature instability. For a better understanding of Eq. (4) one can study the stability of different 
perturbation modes as follows. 

To study the stability of 2D localized structures, we perturb tq with fluctuations of the form ro(0, t) = RoRe[exp{im9+ 
Am^)], where, m is integer number, 6 S [0, 27r[, i?o an arbitrary small constant, and Re denotes the real part. By 
substituting this perturbations in Eq. (4) and linearizing in i?o, we get the growth rate relation 

Kn{rs) = {v- 2/3) . (7) 

Even though the previous analysis has been made considering a fixed value for the radius of the localized structure, 
this radius is determined by balance between the interface energy and the energy difference between the homogeneous 
states. The interface energy and energy difference are proportional to i/ and rj parameter, respectively. The radius 
of the localized structures Tg is proportional to vjri [36], however no analytical expression is available. Hence, the 
possible radius only lower bounded. This lower bound is determined by the size of the core of the front between 
homogeneous states. For the sake simplicity, to figure out the role of the growth rate relation we consider that Vg is 
an independent parameter. Then, for studying the effect of the curvature on the modal instability, one can consider 
an artificial circular structure of any given radius and evaluate its stability. Figure 4 shows the growth rate Am as 
a function of the radius of circular 2D localized structures for different values of angular index m. For fixed values of 
all parameters, the stability of the localized structure is affected with its radius. For small radius, the mode m = 0 
causes radial growth without any change on the structures shape. Only modes m ^ 0 affect the circular shape of the 
localized structures. The angular mode m = 2 becomes unstable (A 2 > 0) and leads to an elliptical deformation of 



FIG. 3. (Color online) Bifurcation diagram of Eq. (1) in (e, rj) space for v = 2.0. In //-zone (dashed black), stable circular 
localized structures are observed. In /-zone (grey crosses) generation of labyrinthine structures are observed from localized 
structures. 
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the circular shape of the localized structure as shown in Fig 2.2. This mode has a larger domain of instability and an 
elliptical deformation occurs for ro = 5.7 as shown in Fig. 4. Triangular deformation corresponding to the angular 
index m = 3 occurs for rg = 7.9. It has been shown that for fixed parameter values, different modes emerge by means 
of the localized structure’s radius increment. Nevertheless, if the radius of the localized structure is fixed, one expects 
that by the variation of a control parameter, lets say one could control the emergence of certain unstable modes 
(cf. Fig. 5). In the above results, we have considered that the radius of the localized structure is large, however, 
this approximation is not always valid. Numerical simulations compared to approximate analytical results of the 
growth rate of localized structures is shown in Fig. 4, numerically, an unstable circular spot of arbitrary radius is 
generated by hand. From this graphs we conclude that the approximate analysis of the C operator considering 
as an independent parameter gives a qualitative description of the observed dynamics. Moreover, the approximated 
growth rates coincide in magnitude to numerical observations, and the transition radius between unstable modes are 
consistent. On the other hand, the stability analysis of £ is only accessible numerically. 


Let us fix the value of the radius to rg = 4.0 and to Vg = 7.0, and let now the diffusion coefficient v be the control 
parameter. For these values, the first angular index m mode to become unstable is to = 2 as shown in Fig. 5. From 
this figure we see that when increasing zz, higher order modes become unstable. Comparison between the two different 
values of the radius shows the higher unstable modes appear at lower range of values for ly when increasing the radius 
^ 0 - 



FIG. 4. (Color online) Lines show approximated theoretical growth rate as function of the localized structure radius for 
TO = 0, 2, 3 using formula (7). Points show numerical results for growth rates of modes: m = 2 (circles), to = 3 (squares), 
and TO = 4 (triangle). With parameters: rj = —0.065; e = 2.45; u = 2-, dx = 0.5; dt = 0.03. The insets show the different 
perturbation modes of a localized structure. 



7 




FIG. 5. (Color online) Theoretical growth rate as function of v using formula (7), for parameters: rj = —0.065; e = 2.45. 
Emergence of instability modes for fixed radii, for 1) Vs = 4.0, and 2) Vs = 7.0. 


III. TRANSVERSAL INSTABILITY OF ROD STRUCTURES AND EMERGENCE OF 

LABYRINTHINE PATTERNS 

The SHE Eq. (1) admits a single stripe like solution [37, 38]. In order to evaluate the threshold over which 
transversal oscillations appear, we perform the stability analysis of a rod-like structure, by a method similar to the 
one performed in Ref. [35]. For this purpose we perturb the single stripe solution as u, = m/(^) + W{x^ JCq) where Uf 
is the single stripe solution and ^ = x — Xo{y,t) the relative position, Xq is the field that accounts for the shape and 
evolution of the finger, and W{x, Xq) << 1 is a non-linear correction of a single stripe. Applying this ansatz in Eq. 
(1) at first order in W and applying the solvability condition [12], the following equation is obtained for the dynamic 
of Xo 


where 


dtXo = -A'dyyXo + ep'd^XoidyXo)^ - d^Xo, 


(d^ufld^uf) ’ 


and A' = (z/ 


2/3'). 


( 8 ) 

(9) 


Thus Xq satisfies a nonlinear diffusion equation. This equation describes the dynamics of an interface between two 
symmetric states [39, 40]. This model is well known for exhibiting a zigzag instability. Analogously, to the previous 
section, when A < 0 the single stripe solution is stable, and for A > 0, the solution is unstable as result of the 
curvature instability. From equation (8) one expects to observe the single stripe becomes unstable by the appearance 
of an undulation. Figure 6 shows the emergence of this undulated rod-like structure. Note that similar dynamical 
behavior is observed in the propagation of cholesteric finger in liquid crystals [15, 16]. Later, this undulated stripe 
is replaced by the emergence of facets that form a zigzag structure. However the higher nonlinear terms control the 







evolution of the single stripe, then the dynamics of initial zigzag is replaced by the growth of undulations without 
saturation as it is depicted in Fig. 2.4. Therefore, the system displays the emergence of a roll-like pattern which is 
formed in the middle section of the structure and invades the system generating invaginated structure, (see Fig. 2.5). 



FIG. 6. (Color online) Transversal instability of a single infinite stripe of Eq. (1). Image shows a section of a 256x256 point 
simulation with periodic boundary conditions using a pseudo-spectral code. Parameters: rj = —0.065, e = 2.45, v = 2, dx = 0.5, 
and dt = 0.03. 


IV. CONCLUSIONS 


In this paper we have described the stability of localized spot in a Swift-Hohenberg equation. We first constructed 
a bifurcation diagram showing different solutions that appear in different regimes of parameters. Then we have shown 
that the angular index m = 2 becomes first unstable as consequence of curvature instability. This instability leads to 
an elliptical deformation of the localized spot. Other type of deformations corresponding to m = 3, 4 and 5 occur also 
in a Swift-Hohenberg equation. We have shown also that for a fixed values of the localized spot radius, the angular 
index m = 2 becomes unstable for small diffusion coefficient while higher order modes m > 2 become unstable for 
large values of a diffusion coefficient. 

When angular index m = 2 becomes unstable, the curvature instability of localized spot produces an elliptical 
deformation leading to the formation of rod-like structure. Subsequently, it generates undulations in the rod-like 
structure. In the course of time, the space time dynamics leads to the formation of invaginated labyrinth structures. 
To understand this dynamics, we have performed the stability of a single stripe localized structure. 

It should be noted that by an offset transformation, u —>■ u + uq, where ug is a constant, Eq. (1) can be rewritten in 
such a way that the constant parameter rj is removed and a quadratic nonlinearity appears. This quadriatic model is 
equivalent to Eq. (1). The model with a quadratic nonlinearity has been well studied (see the textbook [12] and the 
references therein). This equivalence implies that the results of the present work are also valid for physical systems 
described by the quadratic model. 
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